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Abstract 



We report measurements of double radiative capture in pionic hydrogen 
and pionic deuterium. The measurements were performed with the RMC 
spectrometer at the TRIUMF cyclotron by recording photon pairs from 
pion stops in liquid hydrogen and deuterium targets. We obtained abso- 
lute branching ratios of (3.02 ± 0.27(stat.) ± 0.3l(syst.)) x 10~ 5 for hydro- 
gen and (1.42 if};^ (stat.) ± 0.11(syst.)) x 10~ 5 for deuterium, and relative 
branching ratios of double radiative capture to single radiative capture of 
(7.68±0.69(stat.) ±0.79(syst.)) x 10~ 5 for hydrogen and (5.44 ±g|| (stat.) ± 
0A2(syst.)) x 10 -5 for deuterium. For hydrogen, the measured branching ratio 
and photon energy-angle distributions are in fair agreement with a reaction 
mechanism involving the annihilation of the incident ir~ on the 7r + cloud of 
the target proton. For deuterium, the measured branching ratio and energy- 
angle distributions are qualitatively consistent with simple arguments for the 
expected role of the spectator neutron. A comparison between our hydrogen 
and deuterium data and earlier beryllium and carbon data reveals substantial 
changes in the relative branching ratios and the energy-angle distributions and 
is in agreement with the expected evolution of the reaction dynamics from an 
annihilation process in S-state capture to a bremsstrahlung process in P-state 
capture. Lastly, we comment on the relevance of the double radiative process 
to the investigation of the charged pion polarizability and the in-medium pion 
field. 
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I. INTRODUCTION 



A cornerstone of quantum chromodynamics at low energies is the pion's 
emergence as the approximate Goldstone boson of a spontaneously broken 
chiral symmetry. Consequently, the pion is central to understanding the 
realization of QCD symmetries and their associated currents in low energy 
hadronic processes. In particular, the emission and absorption of pions by 
nucleons is closely connected to the low energy manifestation of the partially 
conserved axial current. 

Radiative pion capture, ttN — > jN, and pion photo-production, 7./V — > 
irN, are well-known benchmarks for experimental tests of the low energy re- 
alization of approximate chiral symmetry. The chiral predictions for pion 
photo-production and radiative pion capture were originally obtained from 
current algebra techniques [1] and are nowadays derived from heavy baryon 
chiral perturbation theory [2,3]. Such arguments relate the threshold am- 
plitude for absorption or emission of s-wave pions to the axial currents of 
the nucleon and the pion. The resulting relationship for ^yp — > ir + n photo- 
production - the Kroll-Ruderman theorem and higher-order corrections - is 
experimentally verified to a few-percent accuracy. 

Less well known - both theoretically and experimentally - is double ra- 
diative pion capture, irN — > jjN, and radiative pion photo-production, 
7 AT — > jirN. Nonetheless the absorption or emission of pions in double radia- 
tive capture or radiative photo-production is similarly constrained by chiral 
symmetry. Predictions for double radiative pion capture ir~p — > 7771, have 
been derived from a current algebra approach by Lapidus and Muskhanov [4] 
and a chiral Lagrangian approach by Beder [5,6] and Cammarata [7]. Most 
interestingly, all studies have suggested that threshold absorption from i = 
orbitals involves the annihilation of the incident tt~ on the tt + cloud of the tar- 
get proton, the double radiative process thus exhibiting intriguing sensitivity 
to the nucleon's pion cloud in the small-to-moderate q 2 regime. 

Originally, some excitement was generated by the prediction of it~it + — > 
77 annihilation in double radiative capture. The context was speculations 
on pion condensation or precursor effects in ordinary nuclei, which motivated 
Ericson and Wilkin [8], Nyman and Rho [9], and Barshay [10] to advance the 
(-7T, 27) reaction on nuclear targets as a selective probe of the pion field in the 
nuclear medium. While pion condensation at nuclear densities is nowadays 
unimaginable, the investigation of the pion field in the nuclear medium is still 
of continuing importance. 

Recently, new studies [11] have attempted to determine the pion's elec- 
tromagnetic polarizability via radiative pion photo-production. The method 
involves the Compton scattering of the incident 7-ray from the 7r + cloud of the 
target proton - i.e. the proton acting as a pion target - with the polarizability 
being encoded in the amplitude for the jtt — > jtt scattering. Consequently, 
a better understanding of the pion cloud contribution to 77r — > jit scattering 
in radiative photo-production or tttt — > 77 annihilation in double radiative 
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capture is of some importance. 

Unfortunately, the predicted branching ratios for double radiative cap- 
ture in pionic atoms are small - of order of 10 -4 - and the experimental 
work is extremely challenging. Nevertheless, experiments by Deutsch et al. at 
CERN [12] and Mazzucato et al. at TRIUMF [13] were successful in record- 
ing gamma-ray pairs from double radiative capture on 9 Be and 12 C targets. 
However, the elementary process of double radiative pion capture on hydro- 
gen, and the simplest nuclear process of double radiative pion capture on 
deuterium, have so far remained unobserved. 1 

Herein we report the results of experimental work on double radiative 
capture in pionic hydrogen and pionic deuterium. The experiment was con- 
ducted with the RMC spectrometer at the TRIUMF cyclotron and detected 
photon-pairs from pion stops in liquid hydrogen and deuterium targets. Our 
results include absolute branching ratios, relative branching ratios, and pho- 
ton energy-angle distributions for double radiative capture on hydrogen and 
deuterium. Our goals in conducting the measurements were to: test the pre- 
dictions for the reaction dynamics of the elementary process, compare the 
results for the elementary process and the simplest nuclear process, and ex- 
amine the evolution of double radiative capture from Z = 1 to complex nuclei. 
A letter on our experimental results for the elementary process has been pre- 
viously published [15]. 2 

This article is organized as follows. In Sec. II we briefly outline the the- 
oretical aspects of the irN — > 77./V process. In Sec. Ill we discuss the back- 
ground difficulties arising from real 7-ray coincidences and accidental 7-ray 
coincidences from pion charge exchange, single radiative capture and other 
processes. Sec. IV describes the experimental setup and trigger electronics. 
The identification of the double radiative capture events, subtraction of the 
two-photon background events, and determination of the two-photon detec- 
tion efficiency are described in Sec. V. In Sec. VI we discuss our results for 
hydrogen and deuterium, and compare our experimental data for these ele- 
mentary processes with earlier data for the nuclear process. Our conclusions 
are presented in Sec. VII. 



X A branching ratio upper limit < 5.5x10 4 for double radiative capture on pionic hydrogen was 
obtained by Vasilevsky et al. [14]. 

2 In addition, in Refs. [16,17] we conducted searches for exotic baryon and dibaryon resonances by 
detection of photon-pairs following pion capture in hydrogen and deuterium. 
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II. THEORETICAL OVERVIEW 



A. Initial state 

Pionic atoms are initially formed in atomic orbitals with large principal 
quantum numbers (n) and an approximately statistical population of the or- 
bital sub-states (£). The atoms then de-excite by a combination of Auger 
transitions and radiative transitions with the cascade ultimately terminating 
in the nuclear absorption process (for details see Ref. [18]). 

In pionic hydrogen and pionic deuterium the nuclear absorption is domi- 
nated by S-state capture from n = 3 — 6 orbitals [19-22]. This unique feature 
of Z = 1 capture originates from the combination of strong Stark mixing be- 
tween degenerate ^-orbitals and weak I > nuclear absorption and is firmly 
established by the available pionic x-ray and strong interaction data [19-23]. 

By comparison, in pionic beryllium and pionic carbon the nuclear cap- 
ture is predominantly from P-states. For 7r~ 12 C, Sapp et al. [24] have given 
summed absorption probabilities of 0.08 ±0.03 for S-state capture, 0.92 ±0.03 
for P-state capture, and ~7x 10~ 4 for D-state capture. For 7r" 9 Be, Sapp et al. 
[24] have tabulated x-ray yields that imply summed absorption probabilities 
of ~0.2 for S-state capture and ~ 0.8 for P-state capture. 

B. Reaction kinematics 

For the ir~p — > 7771 elementary process the important kinematical quan- 
tities are shown in Fig. 1. We denote the pion 4-momentum ir^ = (m n , 0), the 
proton 4-momentum p^ = (M p ,0), the neutron 4-momentum = (E n ,p~^), 
and the photon 4-momenta as = (a>i,fci) and k% = (^2,^2)- The pion, 
proton and neutron masses are m n , M p and M n respectively. 

The measured quantities in our experiment are the two photon energies 
cj±, L02 and the photon-pair opening angle cosine, y = cos 9. In terms of 
measured quantities the 3-momentum transfer q = p" n and energy transfer 
q Q = E n - M p are 

\q\ = \Ju\+ujI + 2uJiUJ 2 y 

and 

q Q = m n - ui - UJ2 

where the binding energy of the pionic atom has been neglected. Note that 
the 4-momentum transfer-squared q 2 is space-like and ranges from zero, for 
back-to-back photons with equal energy partition, to -m 2 , for either parallel 
photons or a photon energy at the end-point energy. 

In the ir~p — > 77/1 elementary process the photon energies range from 
to 129.4 MeV and the neutron kinetic energy ranges from to 8.9 MeV. 
The 7771 3-body phase space distribution strongly favors the equal sharing 
of the available energy between the two photons and somewhat favors large 
photon-pair opening angles over small photon-pair opening angles. 



C. Reaction dynamics 



Calculations of the n~p — > 7771 elementary process have been performed 
by Joseph [25], Lapidus and Musakanov [4] and Beder [5] for at-rest capture 
and Beder [6] and Cammarata [7] for in-flight capture. For at-rest capture on 
hydrogen the predicted absolute branching ratios are 5.1 xlO -5 [25], 5.1 xl0~ 5 
[4] and 5.4xl0 -5 [5], and the corresponding ratios ^(27/17) of double radia- 
tive capture to single radiative capture are 1.4xl0~ 4 , 1.4xl0" 4 andl.3xl0" 4 , 
respectively. 

At tree-level the 7^7 effective Lagrangian of Beder [5] implies three cat- 
egories of contributions to double radiative capture: the tttt diagrams of Fig. 
2.a-b, with both photons attached to the pion, the NN diagrams of Fig. 2.c-e, 
with both photons attached to the nucleon, and the irN diagrams of Fig. 2.f-g, 
with one photon attached to the pion and one photon attached to the nucleon. 
The graphs in Fig. 2.a-b and 2.f-g are leading order O(l), and the graphs in 
Fig. 2.c-e are next-to-leading order 0(1/M n ), in the modern terminology of 
heavy-baryon chiral perturbation theory [28]. 

In threshold i = capture {e.g. X H and 2 H) the contributions of dia- 
grams 2.b and 2.f-g vanish and the only surviving leading-order graph is the 
four-point tttt annihilation graph of Fig. 2. a. The four-point tttt annihilation 
diagram yields a squared matrix element with kinematical dependence [5] 

\M.E.\ 2 oc (u\ + ujI + 2uj 1 u 2 y)(l + y 2 ) / '(^(wi + u 2 ) - wiw 2 (l - y)) 2 

that peaks at small opening angles and equal energy partition. This sole 0(1) 
contribution to I = capture, involving annihilation of the stopped 7r~ on 
the 7r + -field of the target proton, and offering sensitivity to the nucleon's pion 
cloud and the tttt — > 77 vertex, is the most intriguing feature of the double 
radiative process. 

In addition, the NN diagrams of 2.c-e introduce a next-to-leading order 
0(1/ M n ) contribution to threshold £ = capture. These graphs yield a 
squared matrix element with kinematical dependence [5] 

\M.E.\ 2 oc (1 -y) 2 

that peaks at large opening angles but is independent of the energy partition 
between the two photons. According to Beder [5], the 0(1) tttt graph 2. a 
contributes roughly 65%, the 0(1/M n ) NN graphs 2.c-e contribute roughly 
20%, and their interference terms contribute roughly 15%, to the I = total 
branching ratio. 

In threshold 1=1 capture (e.g. 9 Be and 12 C) the leading order 7riV dia- 
grams 2.f-g are found to dominate, with small contributions from the O(l) tttt 
graphs and negligible contributions from the 0(1/ M n ) NN graphs. The ttN 
bremsstrahlung graphs 2.f-g yield a squared matrix element with kinematical 
dependence [5] 

\M.E.\ 2 oc 1/u 2 + 1/lj 2 + (1 + y 2 )/(2cjifj 2 ) 



that diverges for small photon energies. A striking difference between I = 
capture and 1 = 1 capture is the prediction of a roughly ten-fold increase in 
the ratio #(27/17). 

In addition to the above contributions from the 7riV"7 effective Lagrangian, 
the possible effects of virtual 7r°'s and intermediate A's have been investigated 
by several authors [4,8,10]. The contribution of virtual 7r°'s, which involves 
virtual it * — ► 77 decay following off-shell n~p — > Tr°*n charge exchange, was 
found by Lapidus and Musakanov [4], Ericson and Wilkin [8] and Barshay 
[10] to be entirely negligible. The contribution of virtual A's, which involve A 
excitation by pion absorption and A decay by photon emission, was estimated 
by Beder [5] to be roughly 5% in the t = process. 

D. Medium effects 

Of course, in double radiative capture on nuclear targets the dynamics of 
the elementary process must be folded with the nuclear response. Although 
no calculation for double radiative capture on deuterium has been performed 
a number of studies of double radiative capture on light nuclei have been 
published. The most detailed treatments of nuclear double radiative capture 
were published by Christillin and Ericson [26] and Gil and Oset [27]. 

Christillin and Ericson [26] combined an effective interaction incorpo- 
rating both 7T7T and 7rN graphs with a phenomenological nuclear resonance 
incorporating the giant dipole and other resonances to calculate 12 C dou- 
ble radiative capture. They found that an important effect of the nuclear 
medium was a substantial quenching of the partial rate in the kinematical 
region (k\ + A^) — > due to Pauli blocking. Gil and Oset [27] combined an 
effective interaction that includes the tttt graphs but omits 7rN graphs with 
the Fermi gas model and the local density approximation to calculate 9 Be 
and 12 C double radiative capture. In contrast to Christillin and Ericson they 
found a significant enhancement of the partial rate in the kinematical region 
(ki + £2) — > 0, and commented that medium modification of the pion field has 
a large effect for these kinematics. Note that both calculations indicated a 
large increase in the ratio 11(27/17) for the nuclear process over the elemen- 
tary process. 

III. EXPERIMENTAL CHALLENGES 

The measured values of the branching ratios for the known modes of pion 
absorption on hydrogen and deuterium are listed in Table I. For hydrogen the 
capture modes are pion charge exchange, ir~p — > 7r°n, single radiative cap- 
ture, ir~p — > jn, and radiative pair production, w~p — > e + e~n [29-31]. For 
deuterium the capture modes are non-radiative capture ir~d — > nn, radiative 
capture ir~d — ► "ynn, and pion charge exchange ir~d — > ir°nn [32,33]. By 
comparison the anticipated branching ratios for double radiative capture on 
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hydrogen and deuterium are of order 10 -4 - the tiny branching ratios, neu- 
tral particle final-states and continuum energy-angle distributions together 
making detection quite difficult. Moreover, these difficulties are magnified 
by intense backgrounds from real coincidences and accidental coincidences of 
photon pairs. 

Neutral pion decay tt° — > 77 following at-rest pion charge exchange yields 
real 7-7 coincidences, the tt° recoil energy of 2.9 MeV in ir~p — > Tr°n and 
— 1 MeV in ir~d — > Tr°nn imparting a Doppler broadening on the labo- 
ratory energy-angle distributions of the 7r° — > 77 photons (the two-photon 
background from in-flight charge exchange was negligible compared to at- 
rest charge exchange). In hydrogen the broadening yields an energy spread 
55.0 < uj < 83.0 MeV and an angular spread — 1 < y < —0.91. In deuterium 
the broadening yields an energy spread 61.6 < to < 73.9 MeV and an angular 
spread —1 < y < —0.97. The ir° — > 77 background therefore prohibits a mea- 
surement of double radiative capture at large opening angles and endangers 
the measurement at small opening angles due to the finite resolution of the 
detector system. Since the branching ratio for charge exchange is much larger 
in hydrogen (0.607 ± 0.004 [29]) than deuterium ((1.45 ± 0.19) x 10~ 4 [32]) 
the 7r° — > 77 background is more worrisome in hydrogen than in deuterium. 
Herein we denote the real 7-7 coincidences from ir° — > 77 decay as the "7r° 
background" . 

Charge exchange and single radiative capture yield photon-pairs via ac- 
cidental 77 coincidences. In our experiment the incident beam had a micro- 
structure with a pulse period of 43 ns and a pulse duration of 2—4 ns. Conse- 
quently, the probability of two pions stops in one beam pulse was not negligi- 
ble. Such multiple pion stops can generate accidental 7-7 coincidences via the 
detection of one photon from one pion stop and another photon from another 
pion stop in the same beam pulse. These accidental 7-7 coincidences may 
comprise a pair of photons that originate from: two charge exchange events, 
two single radiative capture events, or one charge exchange event and one sin- 
gle radiative capture event. Such accidental 7-7 coincidences have summed 
energies 106 < uj < 258 MeV and opening angles — 1 < y < +1. Herein we 
denote the accidental 7-7 coincidences from multiple pion stops in one beam 
pulse as the u 2ir background" . 

Lastly, two possible sources of two-photon backgrounds were (i) Dalitz 
decay following pion charge exchange and (ii) fj,~ decay from in-beam muon 
contamination. In background process (i) the decay ir° — > je + e~ can yield 
a photon-pair by identification of one true photon (the Dalitz-decay 7-ray) 
and one fake photon (the Dalitz-decay e + e~ pair). In background process 
(ii) the pT decay yields single photons via internal bremsstrahlung or exter- 
nal bremsstrahlung and photon pairs by accidental coincidences with other 
photons from pion stops. Although both backgrounds are much rarer than 
the 7T° background or the 2tt background, they are less straightforwardly dif- 
ferentiated from the double radiative capture events by the applied cuts (see 
Sec. VA). Herein we denote background (i) as the Dalitz background and 
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background (ii) as the fiir background. 



IV. EXPERIMENTAL SETUP 

Our measurements of (7r~,27) capture were conducted with the RMC 
detector on the M9A beamline at the TRIUMF cyclotron. The RMC detector 
is a photon-pair spectrometer with large solid angle coverage and moderate 
energy resolution. The detector was developed for studies of radiative muon 
capture on hydrogen and nuclei, and has been discussed elsewhere [34,35]. 

The M9A beamline [36] provided a negative pion flux of typically 7 x 
10 5 s" 1 with a central momentum 81.5 MeV/c, momentum bite 9%, and spot 
size 3.4x4.4 cm 2 . The beam had a microscopic time structure consisting of 
a 2-4 ns pulse width with a 43.3 ns pulse separation and a 99% macroscopic 
duty cycle. For the hydrogen measurement we conducted the experiment 
with the M9A radio frequency separator which yielded a ir/e/n beam compo- 
sition of 0.73/0.18/0.09. For the deuterium measurement we conducted the 
experiment without the separator which yielded a n/e/fi beam composition 
of 0.06/0.89/0.05. 3 

The incident beam particles were counted in a 4-element plastic scintilla- 
tor telescope with a total thickness 0.635 cm and stopped in either a liquid 
hydrogen target or a liquid deuterium target. Both targets had cylindrical 
cross sections with dome-shaped front faces, and were located at the geomet- 
rical center of the detector system. The H2 cell was 15 cm in length, 16 cm 
in diameter, and 2.7 t in volume, with Au front and side walls of thickness 
0.25 mm. The D2 cell was 15 cm in length, 15 cm in diameter, and 2.5 I in 
volume, with Cu front and side walls of thickness 0.50 mm. For the H2 exper- 
iment the hydrogen filling had a deuterium contamination of order 10~ 6 and 
for the D2 experiment the deuterium filling had a hydrogen contamination of 
order 10~ 3 . 

The outgoing photons were detected by pair conversion in a cylindrical Pb 
converter and e + -e~ tracking in two cylindrical tracking chambers. The Pb 
converter had a thickness 1.12±0.01 mm, radius 13 cm, and length 35 cm, and 
was centered at the geometrical center of the detector system. The tracking 
chambers comprised a low mass cylindrical multi-wire proportional chamber 
and a large volume cylindrical drift chamber. The MWPC had an inner radius 
26.5 cm, outer radius 27.5 cm, and active length ~75 cm, and consisted of 
768 axial anode wires between 2x384 helical inner and outer cathode strips. 
It provided a single (x, y, z) coordinate for each track. The drift chamber had 
an inner radius 30.2 cm, outer radius 58.0 cm and active length ~90 cm, and 
comprised four concentric layers of drift cells with six anode wires per drift cell. 
Layers one, two and four were axial and provided (x, y)-coordinate tracking 



3 This running condition resulted from the failure of the RF separator in the D 2 experiment. 



8 



information. Layer three was stereo and provided the ^-coordinate tracking 
information. The entire detector system was located in a large solenoidal 
magnet that provided a 1.2 kG axial magnetic field. 

Concentric rings of segmented plastic scintillators were employed for fast- 
triggering. These scintillators comprised: the A/A'-rings, located just inside 
the Pb converter radius, the C-ring, located just outside the Pb converter 
radius, and the D-ring, located just outside the drift chamber radius. The 
A/A'-rings had four- fold segmentation and inner radii of 12.2 and 12.6 cm 
respectively. They were used to veto any charged particles that emerged from 
the vicinity of the target. The C-ring had twelve-fold segmentation and an 
inner radius 15.8 cm. The C scintillators were used to identify the e + e~ pairs 
from 7-ray conversion in the lead converter. The D-ring had sixteen-fold 
segmentation and an inner radius ~59.4 cm. The D scintillators were used 
to identify the individual electron and positron tracks that emerged from the 
drift chamber. 

The drift chamber times were recorded by LeCroy 1879 pipeline TDCs 
while the MWPC hits were recorded by custom-built electronics [37]. The 
pulse height and timing information from the beam counters and trigger scin- 
tillators were recorded by CAMAC ADCs and TDCs. The data acquisition 
was based on two SLAC scanner processors (SSPs) [38] and a VAX 3000 work- 
station. One SSP handled the event readout from the hardware to an internal 
memory, the other SSP handled the data transfer from the internal memory 
to the computer, while the computer handled the data storage, run control, 
online monitoring, etc. 

A. Trigger system 

The two-photon trigger was based upon the hit multiplicities and hit 
topologies in both the trigger scintillator rings and the drift chamber cells. 
The trigger was designed for efficient identification of photon pairs from dou- 
ble radiative capture and efficient rejection of single photons and coincident 
photons from single radiative capture and pion charge exchange. The trigger 
was implemented in four stages that involved: trigger scintillator multiplic- 
ities (level one), trigger scintillator topologies (level two), drift chamber cell 
multiplicities (level three), and drift chamber cell topologies (level four). 

Level one was based on the multiplicities of the hits in the trigger scin- 
tillators. The level-one condition was A-A'->2C->3D, where A- A' imposed a 
veto on charged particles, >2C required a pair of 7-ray conversions, and >3D 
required a minimum of three tracks to fully traverse the drift chamber. 

Level two was based on the topologies of the hits in the C-D trigger scintil- 
lators rings and was implemented using a LeCroy 2372 memory look-up unit 
(MLU) and a LeCroy 4508 programmable logic unit (PLU). First the MLU, 
using the recorded distribution of C-counters hits, defined the allowed distri- 
bution of D-counter hits that was consistent with the topology of a double 
radiative capture event. Then the PLU, using the allowed distribution of D- 
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counter hits, determined the number of so-called "valid Ds", i.e. the number 
of matches between allowed D-counters and fired D-counters. A total of > 2 
non-adjacent C-counter hits were required to define an allowed distribution 
of D-counter hits, and a total of >3 valid Ds were required to fulfill the level 
two trigger requirement. 

Note that two different versions of the level two trigger were implemented 
for the hydrogen experiment and the deuterium experiment. In the hydrogen 
experiment, if two C-counters were back-to-back, i.e. were located on the 
opposite sides of the C-ring, no allowed distribution of D-counter hits was 
defined. This permitted a reduction in the trigger rate due to photon pairs 
from 7r° — > 77 decay following ir~p — > Tr°n charge exchange. In the deuterium 
experiment, the branching ratio for charge exchange is small and therefore 
rejection of back-to-back C-counters was unnecessary. 

Level three was based on the multiplicities of the cell hits in the drift 
chamber and was derived by custom-built OR boards [39], which indicated if a 
cell was hit, and a custom-built analog hit counter [40], which accumulated the 
distribution of hits. The trigger imposed minimum values on the number of 
the individual cell hits in the second drift chamber layer (denoted 712) and the 
number of contiguous cell clusters in the third and fourth drift chamber layers 
(denoted 713-4). The values of ri2 and 713-4 were chosen for good efficiency 
in accepting photon-pair events and good efficiency in rejecting single photon 
events. In the hydrogen measurement we collected some data (45%) with 
(ri2, 71.3-4) > (3,3) and some data (55%) with (n-2, 713-4) > (4,6). In the 
deuterium experiment we collected all data with (712,713-4) > (4,6). 

Level four was based on the topologies of the cell hits in the first layer 
of the drift chamber. This trigger was implemented in software by the SSP 
and was based on cell information from the OR boards. It rejected any event 
with back-to-back cell hits in drift chamber layer one and was useful in the 
reduction of the event rate from ir° — > 77 decay following w~p — > ir°n charge 
exchange. The level four trigger was enabled for the hydrogen measurement 
but disabled for the deuterium measurement 

A dedicated tt° trigger - which allowed the back-to-back topologies of 
7-ray conversions from tt° decays - was used for measurements of the spec- 
trometer acceptance, resolution function, etc. A dedicated beam trigger - 
which generated events from pre-scaled hits in beam counters - was used for 
the determination of the pion content of the incident beam. 

V. DATA ANALYSIS. 

We report hydrogen and deuterium data that were collected in two four- 
week running periods. The hydrogen experiment recorded a total of 3.10x 10 11 
pion stops and the deuterium experiment recorded a total of 3.07 x 10 11 pion 
stops. 

Note that the running conditions for data taking in hydrogen and deu- 
terium were different in several ways. Firstly a separated beam was used in the 
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hydrogen experiment, while an unseparated beam was used in the deuterium 
experiment (see Sec. IV). For deuterium this resulted in a large electron beam 
rate which therefore required greater attention to counting the ir~ stops. Sec- 
ondly the MWPC was not operable at the large electron beam rate of the 
deuterium experiment. For deuterium this resulted in less z-tracking infor- 
mation and poorer z-tracking resolution. Thirdly, the level two trigger and 
the level four trigger were different in rejecting back-to-back hit topologies 
for the hydrogen experiment but accepting back-to-back hit topologies for 
the deuterium experiment. 

The data analysis involved several steps. In step one we applied cuts to 
identify photon pairs that originated from the target. In step two we applied 
cuts to select the photon pairs from double radiative capture and reject the 
photon pairs from various background sources. In step three we determined 
the contributions of any residual background events in the double radiative 
capture events. Lastly, we determined the number of negative pion stops, and 
obtained the photon pair acceptance, in order to derive the branching ratio 
for double radiative capture. 

A. Identification of photon pairs 

The analysis procedure involved first identifying the various tracks in the 
tracking chambers, then pairing the e + e~ tracks to make individual photons, 
and finally pairing the individual photons to make photon pairs. The criteria 
for identifying "good tracks" and making "good photons" are described below. 

1. Tracking cuts 

Each helical track is defined by a circle fit to the (x, y) coordinates and a 
straight-line fit to the (p, z) coordinates that were obtained from the cham- 
ber hits. The tracks fell into two categories: so-called "non-wrap around" 
tracks where the particle emerges through the drift chamber outer radius, 
and so-called "wrap around" tracks where the particle spirals within the drift 
chamber outer radius. We imposed requirements on the maximum variance 
(S% y ,Sz) and the minimum points (N xy ,N z ) in the (x,y) circle fit and (p, z) 
straight-line fit. In addition, for hydrogen, a requirement was imposed on the 
maximum distance d wc of the particle's trajectory from the associated MWPC 
hit. 

The cut parameters for the hydrogen analysis and the deuterium analysis 
are summarized in Table II. Note that the (x,y) resolution is somewhat 
poorer for positrons than electrons - due to the 2.7° Lorentz angle of the drift 
electrons in the 1.2 kG magnetic field - and therefore the values of S xy are 
larger for e + tracks. For hydrogen a cut was also imposed on the distance d wc 
between the fitted track and the corresponding MWPC hit. 
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2. Photon cuts 



The e + e~ tracks from 7 — > e + e~ conversion must converge at the Pb 
converter radius. Therefore in pairing the e + e~ tracks into reconstructed 
photons we imposed a maximum on the x,y separation (d xy ) and the z sep- 
aration (d z ) of the e + e~ tracks at the Pb radius. In addition, we required 
that the reconstructed photon must have originated from the direction of the 
target cell. Therefore a cut was imposed on the magnitude (r c i ose ) and the 
z-component (z c i ose ) of the closest approach of the photon momentum vector 
to the detector's geometrical center 

The values of cuts on d xy , d z , r c i ose and z c i ose for hydrogen and deuterium 
are summarized in Table III. Note that the z-tracking was poorer in the 
deuterium experiment than the hydrogen experiment, due to the inoperable 
MWPC. This is reflected in the different values for d z and z c i ose in hydrogen 
and deuterium. 

3. Photon pair events 

Good photon pairs were defined as events with four good tracks, i.e. four 
tracks passing all tracking cuts in Table II, and two good photons, i.e. two 
photons passing all photon cuts in Table III. From 3.10 x 10 11 pion stops 
in hydrogen a total of 2.3 x 10 6 photon pairs survived the combination of 
tracking cuts and photon cuts. From 3.07 x 10 11 pion stops in deuterium a 
total of 2.3 x 10 5 photon pairs survived the combination of tracking cuts and 
photon cuts. 

The summed energy {u\ + U2) spectra, individual energy [u) spectra, 
and opening angle cosine (y) spectra, for reconstructed photon pairs from 
hydrogen and deuterium, are shown in Fig. 3. The photon energies u> were 
obtained via the sum of the e + e~ momentum vectors at the converter radius, 
i.e. p 7 = p e + +p e -, and the opening angle y was obtained via the scalar 
product of the two 7-ray momenta, i.e. y = p 7 i • p 7 2 / \p-yi | |p 7 2 1 • Note that 
the e + ,e~ transverse momenta pt were derived from the track radius r in the 
circle fit via pt = erB, and the e + e~ total momentum p was derived using 
the track pitch <j) from the straight-line fit via p = pt/ sin <fi. 

The spectra in Fig. 3 are dominated by real photon coincidences from the 
7r° background and accidental photon coincidences from the 2ir background. 
The 7r° background is most readily seen in the opening angle spectra as events 
with y < —0.7 and the 2tt background is most readily seen in the sum energy 
spectra as events with E > 150 MeV. The different numbers of photon pairs 
from 1 H and 2 H - i.e. 2.3 x 10 6 versus 2.3 x 10 5 - reflects the much larger ir° 
background in the 1 H experiment than in the 2 H experiment. 
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B. Identification of double radiative capture events 

Next we describe the cuts we applied on the beam counter pulse height 
and C-counter timing information to remove the 2tt and \m backgrounds and 
the opening angle to remove the ir° background. These cuts are denoted the 
beam counter amplitude cut, the C-counter timing cut and the opening angle 
cut, respectively. 

1. Beam counter amplitude cut 

To distinguish between single pion stops and multiple pion stops we ap- 
plied a cut on the pulse height information from the four beam counters. First, 
the pulse heights from the eight photo-multipliers viewing the four beam scin- 
tillators were normalized and summed to obtain a single beam counter ampli- 
tude. The resulting amplitude spectra from representative 1 H and 2 H runs are 
shown in Fig. 4. The lower energy peak corresponds to events with single ir~ 
stops while the higher energy peak corresponds to events with two 7r~ stops. 
In the X H spectrum the peak widths were dominated by the energy straggling 
of the incoming pions and the finite resolution of the beam counters. In the 
2 H spectrum the peak widths were further increased by additional pile-up 
with beam electrons. 

Also indicated in Fig. 4 are the corresponding beam counter amplitude 
cuts at Chan. 1069 in hydrogen and Chan. 1150 in deuterium. The cut effi- 
ciencies for passing two-photon events from single ir~ stops were determined 
utilizing the photon pairs from ir° — > 77 decay via both dedicated ir° trigger 
runs and standard two-photon trigger runs. We found an efficiency of 100% 
for the hydrogen data and 92% for the deuterium data. 

2. C counter timing cut 

Two categories of accidental 77 coincidences are capable of circumventing 
the beam counter amplitude cut. They are: (i) fin background events involv- 
ing one prompt 7-ray from a ir~ stop and one delayed 7-ray from an earlier fi~ 
stop, and (ii) tttt background events originating from two separate tt~ stops 
in two adjacent beam pulses. Note that the /U7r background was much larger 
in the unseparated-beam deuterium experiment than in the separated-beam 
hydrogen experiment. 

To remove such backgrounds we applied a C counter timing cut. The 
cut was based on the time difference (tci — tc2) between the two C-counters 
that were fired by the two e + e~ pairs in the two-photon event. For hydrogen 
we used a timing window of —4 < (tci — tc2) < +4 ns and for deuterium 
we used a timing window of —5 < (tci — tci) < +7 ns (the tci — tc2 time 
difference spectrum for deuterium photon pairs is shown in Fig. 5). Their 
corresponding efficiencies for passing real 77 coincidences were determined 
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utilizing the photon pairs from tt° — > 77 decay using both dedicated ir° trigger 
runs and standard two-photon trigger runs. We found an efficiency of 99.5% 
for the hydrogen data and 99.0% for the deuterium data. 

3. Opening angle cut 

Finally, we applied a cut on the opening angle cosine y to remove the 
7r° background (the distributions in opening angle for all photon-pairs from 
hydrogen and deuterium are shown in Fig. 3). In our final analyzes we used 
values of y > —0.1 for hydrogen and deuterium although studies were made of 
the sensitivity of our results to the cut value. Note that (i) the ir° background 
was much larger in the hydrogen measurement and (ii) the opening angle 
resolution was significantly poorer in the deuterium measurement. 

The efficiency of the opening cut for passing events from double radiative 
capture is obviously dependent on the angular distribution for the (it, 77) 
events. The cut's effect on (ir, 77) events is discussed in detail in Sec. VI. 

4- Double radiative capture events 

The resulting distributions of events in individual energy, summed energy 
and opening angle - after the tracking, photon, beam counter amplitude, C 
counter timing and opening angle cuts - are shown in Fig. 6 for hydrogen 
and deuterium. For hydrogen the spectra contain a total of 665 two-photon 
events with 597 events having sum energies <150 MeV and 566 events having 
sum energies 80-150 MeV (the sum energy cut <150 MeV removes the vast 
majority of the tttt background events and the sum energy cut >80 MeV 
removes the vast majority of any Dalitz background events). For deuterium 
the spectra contain a total of 521 two-photon events with 335 events having 
sum energies <150 MeV and 327 events having sum energies 80-150 MeV. 

The summed photon energy spectrum shows a peak at E ~ that 
corresponds to the double radiative capture events. 4 The photon pairs with 
summed energies > 150 MeV are dominated by accidental 77 coincidences 
from unrejected tttt background. Note that the energy distribution of the 
accidental 77 coincidences are different in hydrogen and deuterium - the 
former being dominated by gamma-rays following (7r~,7r°) charge exchange, 
the latter being dominated by gamma-rays following (tt~ , 7) radiative capture. 
The individual photon energy spectrum shows a broad continuum with a low 
energy cut-off at about 20 MeV from the acceptance of the spectrometer. At 
individual photon energies > 100 MeV the spectra are also dominated by the 
accidental 77 coincidences from the tttt background. 



4 The peak centroid is shifted downwards due to energy loss of the e + e pairs in traversing 
lead converter, trigger scintillators, etc. 
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A typical double radiative capture event is shown in Fig. 7. It has zero 
hits in the A-counter ring, two hits in the C-counter ring, and four hits in 
the D-counter ring. The topology of trigger scintillator hits and drift cell hits 
is consistent with diverging e + e _ pairs, and therefore the event fulfilled the 
two-photon trigger. 

C. Subtraction of background events 

A small fraction of tt° background events and tttt background events were 
found to survive the beam amplitude, C-counter timing and opening angle 
cuts - and therefore to contaminate the sum energy region <150 MeV of 
double radiative capture events. Below we describe the methods we used to 
subtract these backgrounds. 

1. tt° background 

To estimate the number of tt° background events that survived the cuts 
we utilized (i) the measured opening angle distribution for tt° events at angles 
y < —0.45 and (ii) the simulated opening angle spectrum for tt° events at all 
angles. 

The Monte Carlo simulation was performed using the CERN GEANT3 
package [41] (see Sec. VE for details). Recall, the true angular distributions of 
7-7 coincidences from ir° decay following at-rest charge exchange on hydrogen 
are limited to —1.0 < y < —0.91 (the 7r°'s originating from charge exchange 
on hydrogen dominate both the hydrogen data-set and the deuterium data- 
set). However, the measured angular distribution of 7-7 coincidences from 
7r° decay shows a significant opening angle tail from the pair spectrometer 
response function. 

Fig. 8 compares the measured and simulated opening angle spectra for 
the region —1.0 < y < —0.6 using a dedicated ir° trigger (Sec. IVA). It 
indicates good agreement between measurement and simulation for ir° events, 
and affirms our confidence in using the simulation to determine the angular 
response of the pair spectrometer and to estimate the ir° contamination in 
the region y > —0.1. 

Fig. 8 also compares the angular distribution from measurement and sim- 
ulation for the region —0.8 < y < +1.0 using the standard 27 trigger (Sec. 
IVA) for the photon-pairs surviving all the applied cuts except the angle cut. 
The simulated data in Fig. 8 imply a small fraction of ir° events with recon- 
structed opening angles y > —0.1 - thus contaminating the region of double 
radiative capture events. To estimate the n° contamination we defined two 
windows: a ir° window of —0.7 < y < —0.45 in which the ir° background is 
dominant, and a (ir, 27) window of —0.1 < y < +1.0 in which the (it, 27) 
signal is dominant. Using these windows the n° contamination with y > —0.1 
was determined via 
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N* X = N* X (N% C /N% C ) 

where N^ X is the number of measured events in the n° window and NJ^ C 
and N^i c are the number of simulated events in the 7r° and (it, 27) windows. 
For hydrogen the procedure yielded a ir° background estimate of 53 ± 30 
events, and for deuterium the procedure yielded a ir° background estimate of 
2 ± 1 events. Note that the errors in the estimate of the ir° background are 
dominated by the statistical uncertainties in the ir° simulation. The simula- 
tion was based on generating a total of 9xl0 7 ir° decays from at-rest charge 
exchange. 



2. tttt background 

To estimate the number of -kit background events that survived the cuts 
we utilized their measured sum energy spectra in Fig. 9. The summed energy 
spectra were obtained by selecting the photon pairs failing the beam counter 
amplitude cut. In the hydrogen spectrum (Fig. 9 solid curve), the lower 
energy peak at ~120 MeV is predominantly two 7-rays from two tt° — > 77 
decays while the higher energy peak at ~180 MeV is predominantly one 7-ray 
from 7r° — > 77 decay and another 7-ray from single radiative capture. In the 
deuterium spectrum (Fig. 9 dashed curve), the single broad peak at ~230 MeV 
originates from two 7-rays from two single radiative capture events. 

Clearly the sum energy of tttt events can exceed the kinematical limit for 
double radiative capture events. We therefore defined an upper threshold E Q , 
above which double radiative capture events are absent but tttt background 
events are present, and computed the number of tttt events in the double 
radiative capture region <150 MeV from number of tttt events in the high 
energy region >E . For hydrogen this method gave 48±8 tttt background 
events below 150 MeV in the double radiative capture sum energy spectrum 
of Fig. 6. For deuterium the same method gave 20±2 tttt background events 
below 150 MeV in the double radiative capture sum energy spectrum of Fig. 
6. The total number of tttt events at all energies were ~100 in hydrogen and 
~200 in deuterium. 

A slight complication in this approach was the finite resolution of the 
pair spectrometer. This generates a small high energy tail of double radiative 
capture events that extends beyond their kinematic limit. We therefore varied 
the upper threshold E Q from 160 to 180 MeV to examine the sensitivity to 
effects of the high energy tail of the spectrometer resolution function. For both 
hydrogen and deuterium the resulting variations in the resulting estimates of 
the tttt contamination were smaller than la. 
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3. Other backgrounds 



One potential source of background events was accidental coincidences 
between a delayed photon from \x decay and a prompt photon from ir cap- 
ture - the [Air background. Such events were present in the unseparated-beam 
deuterium measurement but absent in the separated-beam hydrogen mea- 
surement. In deuterium, the contamination from [iix events that survived the 
C-counter timing cut was estimated from the ^7r continuum background in 
the tci-tc2 time spectrum. This procedure yielded a background estimate of 
(15 ± 4) \m events in the deuterium data set. 

Another potential source of background events was ir° — > r ye + e~ Dalitz 
decay following ir~p — > ir°n charge exchange with the combination of the 
Dalitz pair and the 7 conversion yielding a photon-pair event. Although the 
branching ratio for the 7r° — > je + e~ decay mode is small (1.2%) and the in- 
efficiency of the A, A' veto is tiny (<lxl0 -4 ), it is possible for Dalitz events 
to generate photon pairs with small opening angles. We therefore conducted 
a GEANT3 simulation of Dalitz events which showed that the Dalitz events 
with small opening angles (y > —0.1) had small summed energies (<80 MeV). 
Consequently, a worst-case estimate of the Dalitz background in the sum en- 
ergy range 80-150 MeV was made by attributing all events with <80 MeV, to 
Dalitz background, giving a worst-case contribution of < 2 events in hydrogen 
and < 1 event in deuterium. 

Another potential source of background events is nuclear double radiative 
capture from pion stops in neighboring materials. Using (i) the fraction (15 ± 
1)% of pion stops in materials other than liquid H2/D2 (see Sec. VI), and (ii) 
an estimated branching ratio 5 x 10~ 6 for nuclear double radiative capture 
with E> 25 MeV and y > —0.1 [12,13], we derived upper limits of <10 nuclear 
(it, 27) events in the hydrogen data-set and <20 nuclear (n, 27) events for the 
deuterium data-set. In addition, the fiducial cuts on photon trajectories will 
likely remove a large fraction of such events from pion stops in the target walls, 
beam scintillators, trigger scintillators, etc. We therefore treated the above 
estimates as very conservative upper limits on the nuclear (n, 27) background. 

4- Background- subtracted events. 

A compilation of the various estimates of the background sources in the 
hydrogen experiment and the deuterium experiment are given in Table IV. 
After subtracting the backgrounds, we obtained a total of 465±4q double 
radiative capture events with summed energies from 80 to 150 MeV in hy- 
drogen and a total of 292±3,7 double radiative capture events with summed 
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energies from 80 to 150 MeV in deuterium. 5 The quoted uncertainties include 
the statistical errors and the systematic errors in the signal counts and the 
background estimates. 

D. Counting of incident pions 

The determination of the total number of incident pions was made by 
counting the total number of incident beam particles and measuring the pion 
fraction of the incident beam. The latter measurement was performed uti- 
lizing the beam trigger (Sec. IV A) and resulting amplitude and timing in- 
formation from the beam counters. Note that the beam trigger events were 
collected simultaneously with the 27 trigger events for the duration of the 
measurements. 

Because of the 43 ns micro-time structure of the primary proton beams 
the arrival of secondary pions, muons and electrons are separated in time. At 
81.5 MeV/c the FWHM of pion arrival times was roughly 4.6 ns, while the 
separation from electron and muon arrivals was ~14 and ~21 ns, respectively. 
In addition, the different particles have different energy loss and therefore 
different pulse heights in the beam counters. The combination of arrival time 
and pulse height at the beam counters was sufficient to clearly differentiate 
between incoming particles in both the separated beam and the unseparated 
beam. A representative plot of arrival time versus pulse height for deuterium 
- showing the quality of separation between pions, muons and electrons - is 
shown in Fig. 10. 

Based on the pulse height and time discrimination between particle types 
we determined an average beam composition for the hydrogen experiment of 
73% 7r~, 18% e~ and 9% \i~ and the deuterium experiment of 6% 7r~, 89% e~ 
and 5% pT . These pion fractions were used to derive the number of livetime- 
corrected incident pions in the hydrogen experiment of (3.10±0.03) x 10 11 and 
the deuterium experiment of (3.07±0.15) x 10 11 . The quoted uncertainties are 
based on the analysis of the pion fraction in the incident beam and are larger 
for the unseparated-beam deuterium experiment than for the separated-beam 
hydrogen experiment. 

E. Acceptance of pair spectrometer 

To compare the experimental data to theoretical calculations one requires 
a determination of the two-photon response of the detector system. The 
response function represents the probability of reconstructing photon pairs 



5 Note that in Ref. [15] we quote 482±42 as the number of double radiative capture events with 
summed energies >80 MeV rather than 465±4Q events with summed energies 80—150 MeV. How- 
ever, the data set is identical. 
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with true energies lu±, u)2 and opening angle y with measured energies u^, 
cu' 2 and opening angle y' , thus yielding the measured distribution of photon 
pairs from the true distribution of photon pairs. The response includes the 
detector geometry, the interaction and conversion of photons, the energy loss 
and multiple scattering of electrons, and the effects of the hardware trigger 
and the software cuts. 

To determine the detector response function we employed a Monte Carlo 
simulation (RMCGEANT) that was based on the CERN GEANT3 package 
[41]. It incorporated the full geometry of the detector system - i.e. beam 
counters, target assembly, trigger counters, tracking chambers and the sur- 
rounding magnet - and all relevant interactions of photons and electrons. 
The simulation produced data files containing Monte Carlo events that were 
analyzed with the same software and the same cuts as the measured data. 
The program included event generators for pion charge exchange ir~p — > 7r°n, 
single radiative capture ir~p — > 771 and double radiative capture ir~p — > 7771. 
The generator we used for double radiative capture was based on formulas 
contained in Beder's article [5]. 

An important test of the GEANT simulation of the detector system was 
afforded by the photon pairs from the tttt background. These pairs originate 
from separate ir~p — > 771 and/or ir~p — > TT°n reactions from two pion stops 
in one beam pulse. Since the absolute yields and energy-angle distributions 
of such pairs are well known, they offer an invaluable means of calibration 
across the full range of photon energies and opening angles. In the hydrogen 
experiment this calibration procedure was performed continuously for the 
entire duration of the data taking. In the deuterium experiment the procedure 
was limited to dedicated calibration runs with liquid H2 filling. 6 

Representative measured and simulated spectra for photon-pairs from tttt 
events are compared in Fig. 11. The energy-angle distributions from experi- 
ment and simulation are in agreement at the level of a few percent and affirms 
the high quality of the GEANT simulation of the detector system. Note that 
the absolute efficiency for detecting photon pairs was found to be smaller in 
the experiment than the simulation by a factor F = 0.90 ± 0.09 in hydrogen 
and F = 0.96±0.06 in deuterium. This difference is believed to originate from 
detector inefficiencies that are present in the experiment but are absent from 
the simulation. These inefficiencies include trigger scintillator inefficiencies, 
tracking chamber inefficiencies, and possible effects from chamber noise on 
track recognition. For both hydrogen and deuterium the run-by-run variation 
of the normalization factor was ±4% and the total change over the running 
period was <10%. 



During running on deuterium the random 7-7 coincidences from the tttt background and true 
7-7 coincidences from the tt° background were used to monitor the relative acceptance for photon 
pairs. An absolute determination of the acceptance was not possible due to uncertainties in the 
hydrogen contamination in the liquid deuterium. 
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Another important test of the overall quality of the GEANT simulation 
of the detector system was afforded by photon pairs from ir° — > 77 decay 
following w~p — > 7r°n charge exchange. These data were collected in runs using 
a dedicated 7r° trigger in both the hydrogen experiment and the deuterium 
experiment. Likewise the measured data and simulated data for photon pairs 
from ir° — > 77 decay show good agreement and this also affirms the high 
quality of the GEANT simulation of the detection system. The normalization 
factors of F = 0.90 in hydrogen and F = 0.96 in deuterium that were derived 
from the tttt study are also consistent with the ir° study. 

The detector acceptance eAO for double radiative capture was obtained 
using the RMCGEANT package and the double radiative capture event gen- 
erator. The procedure yielded acceptances of eQ = 0.66 x 10~ 4 for hydrogen 
and eSl = 0.90 x 10 -4 for deuterium when using the standard energy-angle 
cuts and assuming photon-pairs with energy-angle distributions derived by 
Beder [5]. Note that the acceptance for deuterium is larger than hydrogen 
due to the looser cuts that were applied in the presence of the smaller tt° 
background. 



VI. RESULTS 



Herein we discuss the experimental results from the hydrogen experiment 
and the deuterium experiment. We first describe the determination of the 
branching ratios in the two experiments and then discuss the comparison 
between theory and experiment for the elementary process and the nuclear 
process. 



A. Branching ratios for hydrogen and deuterium 

The absolute branching ratios for double radiative capture on pionic hy- 
drogen and pionic deuterium were obtained via 



N tt- ■ • F ■ c s top ■ c brn 

where iV 77 is the number of background-subtracted double radiative capture 
events, N n - is the number of livetime-corrected incident pions, eA$7 • F is the 
detector acceptance, and c s t op and c& m are minor correction factors which are 
discussed below. 

The values employed in computing the branching ratios for hydrogen and 
deuterium are summarized in Table V. For the double radiative capture events 
with sum energies 80 < uj\ + u>2 < 150 MeV and opening angles y > —0.1 we 
took A 77 = 465±4Q for hydrogen and iV 77 = 292±<j,7 f° r deuterium from Sec. 
VC4. For the livetime-corrected pion stops we took = (3.10±0.03) x 10 11 
for hydrogen and = (3.07 ± 0.15) x 10 11 for deuterium from Sec. VD. 
The GEANT simulation of detector acceptances of Sec. VE gave eAQ = 
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0.66 x 10 -4 for hydrogen and eAfi = 0.90 x 10~ 4 for deuterium. The factors 
F = 0.90 ± 0.09 for hydrogen and F = 0.96 =L 0.06 for deuterium account for 
small differences between the measured and the calculated acceptances due to 
various detector efficiencies, etc. The factor c s t op = 0.85 ± 0.01 accounts for 
the fraction of incident pions stopping in liquid hydrogen or liquid deuterium 
and was computed vici ci Monte Carlo simulation [34] . The factor Cb m — 0.98 
for hydrogen and q> m = 0.91 for deuterium accounts for efficiencies of events 
surviving the beam telescope amplitude and C-counter timing cuts that are 
omitted from the simulation (the uncertainties in q, m are negligible). 

Using these values and Eqn. 7 we obtained absolute branching ratios of 
(3.02 ± 0.27(stat.) ± 0.31(syst.)) x 10~ 5 for hydrogen and (1.42 ±g;?| (stat.) ± 
0.11(syst.)) x 10 -5 for deuterium. The corresponding values of ratios #(17/27) 
of double radiative capture to single radiative capture are (7.68=L0.69(stat.) ± 
0.79(syst.)) x 10~ 5 for hydrogen and (5.44±g;| 4 (stat.) ±0.42(syst.)) x 10~ 5 for 
deuterium. The quoted uncertainties are dominated by the statistical error 
in the number of photon pairs and the systematic error in the acceptance for 
double radiative capture events. 

For completeness we also quote the partial widths T for double radia- 
tive capture. We use the recent measurements of the total widths of the IS 
atomic states of 0.86±0.07 eV in pionic hydrogen (Schroeder et al. [42]) and 
1.13±0.13 eV in pionic deuterium (Hauser et al. [43]). Combined with our 
absolute branching ratios they yield partial widths of T±s = 26.0 ±4.1 //eV 
for hydrogen and T±s = 16.2 ± 2.5 /ieV for deuterium. 

The sensitivity of the branching ratios to the opening angle cut and the 
sum energy cut are given in Table VI. Variations of the opening angle cut 
from —0.2 > y > 0.0 changed the branching ratio by 7% in hydrogen and 2% 
in deuterium. Variation of sum energy minimum from 60 to 80 MeV and sum 
energy maximum from 150 to 170 MeV changed the branching ratios by 5% 
in hydrogen and 7% in deuterium. A separate analysis of the two hydrogen 
data-sets with the different level-three triggers (Sec. IV A) gave branching 
ratios that differed by 7%. In all cases the variations are significantly smaller 
than the quoted errors on the branching ratios. 

We stress that the quoted branching ratios are for photon pairs with all 
possible angles and all possible energies. Their derivation assumes the energy- 
angle distributions of photon pairs calculated by Beder [5] - the calculation 
and measurement being consistent for the visible region of the angle spectra 
y > —0.1 and energy spectra to > 25 MeV. 

To test the sensitivity to the energy-angle distribution of the photon pairs 
we also extracted branching ratios with two other energy-angle distributions. 
Assuming an energy-angle distribution given by only tttt annihilation graphs 
the resulting branching ratios are decreased by 13% in hydrogen and 14% in 
deuterium, the decrease reflecting the larger proportion of photon pairs with 
small opening angles. Assuming an energy-angle distribution given by three- 
body 77n phase space the resulting branching ratios are increased by 5% in 
hydrogen and 6% in deuterium, the increase reflecting the smaller proportion 
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of photon pairs with small opening angles. 

B. Discussion of results for hydrogen 

Our measured absolute branching ratio (3.02±0.27(stat.) ±0.31(syst.)) x 
10 -5 for double radiative capture on hydrogen is in fair agreement with 
the predicted values of 5.1xl0 -5 (Joseph [25]), 5.1xl0~ 5 (Lapidus and 
Musakhanov [4]) and 5.4xl0~ 5 (Beder [5]) that embrace the dominance of 
the 7T7T — > 77 annihilation mechanism. The measured branching ratio is defi- 
nitely grossly under-estimated by either an t = NN mechanism or an £ = 
7rN mechanism and definitely grossly over-estimated by an i = 1 7rN mech- 
anism (see Sec. II C and Ref. [5] for details). Intriguingly, the measurement 
is very close to the prediction 3.5xl0 -5 of the O(l) tttt mechanism in the 
absence of the 0(1/M n ) NN correction (we offer no rationale for the omission 
of the 0{1/M n ) corrections from the NN graphs). 

A comparison between the theoretical predictions and the experimental 
data for the two-photon opening angle distribution is shown in Fig. 12. Indi- 
cated are the results of the full i = calculation, O(l) tttt contribution only 
and 0{1/M n ) NN contribution only. Note that the theoretical curves have 
been convoluted with the response function of the detector system and incor- 
porate such effects as energy loss and multiple scattering (see Sec. V E for 
details). The measured spectra for opening angles y > —0.1 are reasonably 
consistent with the angular dependences of both the full calculation and the 
0(1) tttt term only. The experimental data are not consistent with the angu- 
lar dependence of the NN mechanism which falls too steeply with increasing 

y- 

We consider the experimental results for the absolute branching ratio and 
the opening angle distribution to be in support of the underlying theoreti- 
cal idea of a dominant tttt — ► 77 mechanism in the i = elementary process. 
However, the measured branching ratio of 3.02±0.27(stat.)±0.31(syst.) x 10 -5 
and calculated branching ratio of 5.1-5.4xl0 -5 do differ by roughly 40%. For 
radiative pion photoproduction the complete set of leading order (LO), next- 
to-leading order (NLO) and next-to-next-to-leading order (NNLO) graphs 
have been enumerated by Kao et al. [28] within heavy-baryon chiral pertur- 
bation theory; however, no such program has been performed for the double 
radiative process. It is therefore possible that the remaining discrepancy be- 
tween the measured branching ratio and the calculated branching ratio may 
reflect the neglect of either certain NLO graphs or all NNLO graphs in the 
existing calculations of the double radiative process. A modern calculation 
of double radiative capture within heavy baryon chiral perturbation theory 
would be welcome. 

As mentioned in Sec. I, at appropriate kinematics the tttt — > 77 anni- 
hilation graph in double radiative capture and tt^^tt^ scattering graph in 
radiative pion photoproduction are potentially sensitive to the charged pion 
polarizability. Indeed, a recent experiment at MAMI [11], and an earlier ex- 
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periment at Lebedev [44], have attempted to extract the polar izability from 
the irj— >7T7 contribution to the 7p — > 77r + n reaction. Although, for threshold 
double radiative capture the pion polarizability effects are too small to be 
measurable, the dominant contribution of the mr graph and vanishing con- 
tributions of the 7rN graphs do make the process an interesting yardstick for 
related theoretical studies of radiative pion production. 

C. Discussion of results for deuterium 

In comparing our results for hydrogen and deuterium we focus on the ratio 
#(27/17) and the two photon kinematical distributions. A quantitative com- 
parison of absolute branching ratios is more difficult and less meaningful as 
pion absorption on nucleon pairs is important for deuterium but inapplicable 
for hydrogen. 

Using our measured branching ratios for double radiative capture and 
the published branching ratios for single radiative capture we derived values 
for #(27/17) of (7.68 ± 0.69(stat.) ± 0.79(syst.)) x 10" 5 for hydrogen and 
(5.44 ±g;^ (stat.) ±0.42(syst.)) x 1(T 5 for deuterium (Sec. VIA). The results 
imply a decrease in #(27/17) from hydrogen to deuterium by about (29±12)% 
(one small caveat is our detection of photon pairs is limited to y > —0.1 and 
co > 20 MeV). In addition, while the angular distributions are similar for 
hydrogen and deuterium, the energy distributions show evidence of lower 
average photon energies in deuterium (see Fig. 6). 7 

The decrease in #(27/17) and the change in photon energies from hydro- 
gen to deuterium are probably not surprising. Single radiative capture and 
double radiative capture involve different reaction mechanisms and different 
momentum transfer, and therefore effects such as Pauli blocking, initial-state 
wavefunctions, final-state interactions, etc., may be different for these two pro- 
cesses. For example, the lower photon energies in deuterium capture compared 
to hydrogen capture may perhaps reflect the fact that energy is imparted to 
the spectator neutron through the final-state interactions in the deuterium 
process. In addition, the smaller #(27/17) ratio in deuterium capture com- 
pared to hydrogen capture may perhaps reflect a greater effect from Pauli 
blocking at the smaller momentum transfer of the double radiative process. 
Of course, to properly address the hydrogen-deuterium comparison a detailed 
calculation of deuterium capture is needed. 



7 Note that although the photon detection efficiencies were significantly different in hydrogen and 
deuterium their energy dependences were very similar (see Sec. V E) for details) 
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D. Comparison to nuclear double radiative capture. 



The world data for double radiative capture on complex nuclei consist of 
results for 9 Be and 12 C by Deutsch et al. [12] at CERN and results for 12 C from 
Mazzucato et al. [13] at TRIUMF. The experiment of Deutsch et al. obtained 
a partial branching ratio of (1.0 ± 0.1) x 10~ 5 for 9 Be and (1.4 ± 0.2) x 10~ 5 
for 12 C for photon energies exceeding 25 MeV. The experiment of Mazzucato 
et al. obtained a partial branching ratio of (1.2 ± 0.2) x 10~ 5 for 12 C for 
photon energies exceeding 17 MeV. Both experiments reported an increased 
yield of photon pairs at large opening angles although Mazzucato et al. have 
suggested that target contamination by hydrogenous materials may perhaps 
be responsible. The Mazzucato et al. experiment also reported a preference 
for unequal partition of photon energies. 

In comparing the nuclear (tt, 2j) data to our (it, 27) data for hydrogen and 
deuterium we consider the ratio #(17/27) in order to avoid the complication of 
accounting for non-radiative processes in complex nuclei. For 12 C, we derived 
values for #(27/17) of (7.2±1.3) x 10~ 4 from Mazzucato et al. and (8.4±1.3) X 
10~ 4 from Deutsch et al. when normalized using the branching ratio for single 
radiative capture from Bistirlich et al. [45]. For 9 Be, we obtained a value for 
#(27/17) of (4.8±1.3) x 10~ 4 when normalized using the mean branching ratio 
for single radiative capture from the neighboring 7 Li- 10 B nuclei [46]. Clearly - 
as summarized in Table VII - the values of #(27/17) on light nuclei are much 
larger than the corresponding ratios for hydrogen (7.68±0.69±0.79) x 10~ 5 and 
deuterium (5.44 ±g;|| (stat.) ± 0.42(syst.) x 10~ 5 . This increase in #(27/17) 
from Z = 1 to light nuclei was anticipated by several authors [5,26] as an 
evolution from predominantly annihilation dynamics in S-state capture to 
predominantly bremsstrahlung dynamics in atomic P states (see Sec. II C for 
details). For instance, an increase in #(27/17) from £ = capture to t = 1 
capture of approximately ten- fold was predicted by Beder in Ref. [5]. 

Another difference between our data for Z = 1 and earlier data for com- 
plex nuclei is the photon energy spectra. For Z = 1 the energy spectra (Fig. 
6) indicate a preference for the equal partition of the photon energies. For 12 C 
the energy spectra (Rcfs. [12] and [13]) indicate a preference for the unequal 
partition of the photon energies. This difference between t = capture and 
t = 1 capture was also anticipated in early theoretical studies of double radia- 
tive capture by Beder [5] and Christillin and Ericson [26] , and likewise reflects 
the evolving dynamics from annihilation in S-state capture to bremsstrahlung 
in P-state capture. 

Another observation of experiments on nuclei was a backward peak in the 
angular distribution of the photon pairs. This feature was not predicted by 
the earlier calculation of Christillin and Ericson [26] but was obtained in the 
later calculation of Gil and Oset [27]. Unfortunately, the intense background 
of back-to-back photon pairs from tt° — > 77 decay made unthinkable the 
measurement of Z = 1 double radiative capture at large opening angles. 

In summary, we claim the large difference in #(27/17) between the ele- 
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mentary and nuclear processes, and significant difference in the energy par- 
tition between the elementary and nuclear processes, are clear evidence of 
the evolving dynamics of the double radiative process from S-state to P-state 
capture. This conclusion implies that earlier arguments for the sensitivity of 
the (tt~,2j) nuclear reaction to the pion field in the nuclear medium, that 
were founded on the dominance of the tttt — > 77 mechanism, are likely not 
justified. 

VII. CONCLUSION 

We have performed the first measurements of double radiative capture on 
pionic hydrogen and pionic deuterium. The measurements were performed 
using the RMC spectrometer at the TRIUMF cyclotron by recording photon 
pairs from pion stops in liquid hydrogen and deuterium targets. We have 
determined the absolute branching ratios for double radiative capture, the 
ratios #(27/17) of double radiative capture to single radiative capture, and 
the energy-angle distributions of the resulting photon pairs. 

For hydrogen we obtained an absolute branching ratio of (3.02 ± 
0.27(stat.) ± 0.31(syst.)) x 1(T 5 and ratio #(27/17) of (7.68 ± 0.69(stat) ± 
0.79(syst)) x 10~ 5 . The measured branching ratios and energy-angle distri- 
butions are in fair agreement with the theoretical predictions of Joseph [25], 
Lapidus and Musakanov [4] and Beder [5] and lend support to the theoreti- 
cal claim of a tttt^jj annihilation mechanism. This mechanism implies that 
double radiative capture from £ = orbitals inherits an intriguing sensitivity 
to the nucleon's pion cloud and the tttt — > 77 vertex. 

For hydrogen the measured and calculated branching ratio do, however, 
differ by approximately 40%. One possible explanation is the neglect of ei- 
ther certain NLO graphs or all NNLO graphs in the existing calculations of 
the double radiative process. Given the relationship between double radiative 
pion capture ir~p — > 7771, and radiative pion photoproduction 7^ — > 771"+ n 
- and their relevance to the determination of the pion's polarizability - we 
believe that resolving this discrepancy is of some importance. A modern 
calculation of double radiative pion capture within heavy-baryon chiral per- 
turbation theory would be welcome. 

For deuterium we obtained an absolute branching ratio of (1.42 ±012 
(stat.) ± 0.11(syst.)) x 10" 5 and ratio #(27/17) of (5.44 ±g;^ (stat) ± 
0A2(syst)) x 10~ 5 . We found that the deuterium ratio #(27/17) is smaller 
than the hydrogen ratio #(27/17) by (29 ± 12)%. We also found the pho- 
ton energies are shifted somewhat to lower energies in the deuterium data 
compared to the hydrogen data. The decrease in #(27/17) and shift in pho- 
ton energies are possibly revealing the spectator neutron's role via final-state 
interactions and Pauli blocking effects. 

Lastly, we compared our results for hydrogen and deuterium with the ear- 
lier results for beryllium and carbon. The relative branching ratios #(27/17) 
for light nuclei are nearly ten times larger than the corresponding values for 
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the hydrogen isotopes. In addition, the energy distributions are markedly 
different for the light nuclei and the hydrogen isotopes; the former favor- 
ing unequal energy partition with the latter favoring equal energy partition. 
These observations are consistent with the dynamics evolving from a tttt — > 77 
annihilation mechanism for £ = capture to a 7rN bremsstrahlung mechanism 
for £ = 1 capture. This implies that nuclear double radiative pion capture 
may not inherit the high sensitivity to the pion field in the nuclear medium 
as was originally proposed. 
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TABLES 



TABLE I. Compilation of previously published branching ratios for negative pion capture on 
hydrogen and deuterium. 



capture mode 



expt. 



branching ratio 



7T p - 

ir~p 

7T~p —> 
TT~d - 

ir~d - 



n 



nn 



'ynn 



Spuller et al. [29] 
Spuller et al. [29] 
Samios [30] 
Highland et al. [32] 
Highland et al. [32] 
MacDonald et al. [33] 



0.607 ± 0.004 
0.393 ± 0.003 

(2.72 ±0.19) xl0~ 3 
0.739 ±0.010 
0.261 ±0.004 

(1.45 ±0.19) xl0~ 4 
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TABLE II. Summary of tracking cuts for positrons and electrons and "nonwrap arounds" and 
"wrap arounds". The upper rows are the hydrogen values and the lower rows are the deuterium 
values. 



track 
parameter 



nonwrap 



wrap 



nonwrap 



wrap 



hydrogen cuts 



S 2 xy {cm 2 ) 
S 2 (cm 2 ) 

N xy 

N z 
d wc (cm) 



< 0.0048 

< 0.2 

> 9 

> 3 

< 1.0 



< 0.0270 

< 3.0 

> 9 

> 3 

< 1.0 



< 0.0038 

< 0.2 

> 9 

> 3 

< 1.0 



< 0.0210 

< 3.0 

> 9 

> 3 

< 1.0 



deuterium cuts 



Sl y {cm 2 ) 
S 2 (cm 2 ) 



xy 

N z 



d wc (cm) 



< 0.0048 
< 0.2 

> 9 

> 3 



< 0.0270 
< 3.0 

> 9 

> 3 



< 0.0038 
< 0.2 

> 9 

> 3 



< 0.0210 
< 3.0 

> 9 

> 3 
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TABLE III. Summary of the photon cuts for "non-wrap around" photons (i.e. containing no 
wrap-around tracks) and "wrap around" photons (i.e. containing wrap-around tracks). The upper 
rows are the hydrogen values and the lower rows are the deuterium values. 



photon nonwrap wrap 

parameter 



hydrogen cuts 

d xy (cm) < 6.0 < 7.0 

d z (cm) < 6.0 < 7.0 

r c iose(cm) < 8 < 8 

z c iose(cm) -12 < z < +14 -12 < z < +14 

deuterium cuts 

d xy (cm) < 6.0 < 7.0 

d z (cm) < 12.0 < 14.0 

r c iose(cm) < 8 < 8 

z c i ose (cm) -25 < z < +27 -25 < z < +27 
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TABLE IV. Summary of background estimates and background limits for the hydrogen anal- 
ysis and the deuterium analysis. The first row is the two-photon events with summed energies 
80-150 MeV before the background subtraction and the last low is two-photon events with summed 
energies 80-150 MeV after the background subtraction. 



11 1 

background source 


T T 

H 2 


I A 

D2 


80-150 MeV events before sub. 


566 


327 


7T° bkd. 


53±30 


2±1 


7T7T bkd. 


40±8 


20±2 


/j,tt bkd. 





15±4 


Dalitz bkd. 


<2 


<1 


nucl. (it, 27) 


<10 


<20 



80-150 MeV events after sub. 



465±^ 



292±$ 



TABLE V. Summary of quantities for the calculations of the total branching ratios of double 
radiative capture on pionic hydrogen and pionic deuterium. Note that the quoted acceptances 
assumes the energy-angle distributions of Beder. 



quantity 



H 2 



D 2 



eAtt 
F 

Cbm 

absolute B.R. 
ratio R(27/l7) 



465±39 
(3.10±0.03)xlO n 
0.66xl0~ 4 
0.90±0.09 
0.85±0.01 
0.98 

3.02±0.27(stat)±0.31(syst) x 10~ 5 
7.68±0.69(stat)±0.79(syst) x 10~ 5 



292±$ 
(3.07±0.15)xlO n 
0.90xl0~ 4 
0.96±0.06 
0.85±0.01 
0.91 

1.42±[! : ^(stat)±0.11(syst)xl0- 5 
5.44±[j : ^(stat)±0.42(syst)xl0- 5 
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TABLE VI. Sensitivity of the absolute branching ratio to the analysis parameters, different 
data-sets and the theoretical energy-angle distributions for the hydrogen analysis and the deuterium 
analysis. We have combined in quadrature the various statistical and systematic uncertainties. 



cuts / data-set / theory 



H 2 



D 2 



80 < wi +0J2 < 150, y > 0.0 
80 < wi + oj 2 < 150, y > -0.1 
80 < wi + lo 2 < 150, y > -0.2 

70 < uji + uj 2 < 150, y > -0.1 
60 < ui + uj 2 < 150, y > -0.1 

80 < ui + 0J2 < 160, y > -0.1 
80 < ui + 0J2 < 170, y > -0.1 

hydrogen data-set I 
hydrogen data-set II 



2 


94 


± 





41) 


X 


10 


3 


02 


± 





41) 


X 


10 


3 


14 


± 





41) 


X 


10 



;3. 11 ±0.41) x 10~ 5 
;3.17±0.41) x 10" 5 

3.11 ±0.41) x 10" 5 
3.05 ±0.41) x 10" 5 

2.94 ±0.49) x 10~ 5 
3.15 ±0.49) x 10- 5 



Beder, all graphs 


(3.02 ±0.41) 


x 10" 5 


(1.42±8;it) x 


10 


Beder, tttt graphs 


(2.63 ±0.41) 


x 10" 5 


(1.22±°;£) x 


10 


phase space 


(3.16 ±0.41) 


x 10" 5 


n/a 





(i.39±g:iS) x io- 5 

(1.42±g;iS) x 10- 5 
(1.39±g:iS) x io- 5 

(1.48±g;i|) x IO" 5 
(1.48±8:it) x 10- 5 

(i.42±8:it) x io- 5 
(i.38±8it) x io- 5 

n/a 
n/a 



TABLE VII. Summary of determinations of ratios R{2 r )/l r y) from our hydrogen and deuterium 
data and earlier beryllium and carbon data. Note that for beryllium and carbon the ratios are for 
photon energies >25 MeV in Deutsch et al. and >7 MeV Mazzucato et al. 



target 



Ref. 



#(27/l7) 



2 H 

9 Be 
12 C 

12 C 



this work 

this work 
Deutsch et al. [12] 
Deutsch et al. [12] 
Mazzucato et al. [13] 



(7.68 ± 0.69 ± 0.79) x 10~ 5 
(5.44 ± ±8;||(stat.) ± 0.42(syst.) x 10~ 5 
(4.8 ± 1.3) x 10" 4 
(8.4 ± 1.3) x 10" 4 
(7.2 ±1.3) x IO" 4 



33 



FIGURES 




FIG. 1. Definitions of kinematical quantities in the n p ^ 7771 elementary process. 




FIG. 2. Tree-level contributions to double pion radiative capture in the ttNj effective La- 
grangian approach of Beder. Diagrams a-b) are denoted the tttt graphs, c-e) the NN graphs, and 
f-g) 7rN graphs. 
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FIG. 3. The opening angle (top), summed energy (center) and individual energy (bottom) 
spectra of the two-photon events passing the tracking cuts and photon cuts in the hydrogen data-set 
(solid histogram) and the deuterium data-set (dashed histogram) . The tt° background is apparent 
as the events with y < —0.7 in the opening angle spectrum and the tttt background is apparent as 
the events with u\ +ll>2 > 150 MeV in the sum energy spectrum. To facilitate their comparison the 
hydrogen spectra and deuterium spectra are normalized to equal numbers of photon-pair events. 
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FIG. 4. The beam amplitude spectra for the hydrogen data-set photon pairs (lefthand plot) and 
deuterium data-set photon pairs (righthand plot). The plots show the single pion stops associated 
with real photon pairs (the lower amplitude peak) and two pion stops associated with accidental 
photon pairs (the higher amplitude peak). The beam counter cuts, at Chan. 1069. in the hydrogen 
data-set and Chan. 1150 in the deuterium data-set, are indicated by vertical bars. 
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FIG. 5. The time difference (ns) between two C-counters hits for photon pairs in deuterium that 
pass all cuts with the exception of the C-counter timing cut. The events outside the coincidence 
window —5 < (tci — tcz) < +7 ns are dominated by /wr and nir baclground. 
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FIG. 6. The opening angle (top), summed energy (center) and individual energy (bottom) 
spectra for photon pairs after all cuts for the hydrogen data-set (solid histogram) and deuterium 
data-set (dashed histogram). The hydrogen spectra contain a total of 665 photon-pair events with 
y > —0.1 The deuterium spectra contain a total of 521 photon-pair events with y > —0.1 
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RUN 8440 
FVFNT 

FIG. 7. A typical double radiative capture event showing the trajectories of the two e + e~ pairs, 
the circle fits to the four e + e~ tracks, and the fired C and D counters. 
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FIG. 8. Comparison of Monte Carlo and experimental data for photon-pairs from ir° decay for 
the 77 trigger and —0.8 < y+ 1.0 (lefthand plot) and the ir° trigger and — 1.0 < y — 0.6 (righthand 
plot). The solid histograms are the Monte Carlo results and the data points are the experimental 
results. (The apparent difference between Monte Carlo and experimental data for y < —0.75 in 
the 77 trigger data is a consequence of a loose y-cut that was applied in the initial processing of 
the experimental data ) 
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FIG. 9. The photon sum energy spectrum for hydrogen (solid curve) and deuterium (dashed 
curve) for events that fail the beam amplitude cut, i.e. accidental 77 coincidences due to multiple 
pions in a single beam pulse. In deuterium the accidental coincidences are dominated by single 
radiative capture. In hydrogen the accidental coincidences originate from pion charge exchange 
and single radiative capture. 
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FIG. 10. Beam counter amplitude versus time-of-flight for incoming beam particles in the deu- 
terium experiment. The pions are top-left, electrons are bottom-center and muons are right-center. 
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FIG. 11. Comparison of simulated data (solid line) and measured data (data points) for tttt 
background events. The plots are the individual photon energies (upper left), summed photon 
energies (upper right), opening angle (lower left) and energy partition (lower right). Also show is 
the difference between the simulated data and the measured data. 
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FIG. 12. Comparison of the opening angle distribution from the experimental data (closed 
circles) and the theoretical calculations (curves) for double radiative capture in pionic hydrogen. 
The solid curve is the full calculation, the dashed curve is the irir term only, and the dotted curve 
is the NN term only. Note the theoretical curves have been normalized to the summed counts of 
the experimental spectrum for the comparison of the angular dependences. 
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FIG. 13. Comparison of the opening angle distribution from the experimental data (closed 
circles) and the theoretical calculations (curves) for double radiative capture in pionic deuterium. 
The solid curve is the full calculation, the dashed curve is the mr term only, and the dotted curve 
is the NN term only. Note the theoretical curves have been normalized to the summed counts of 
the experimental spectrum for the comparison of the angular dependences. 
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